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Binary snow maps and fractional snow cover data are provided routinely from MODIS (Moderate Resolution 
Imaging Spectro radiometer). This paper investigates how the wide observation angles of MODIS influence 
the current snow mapping algorithm in forested areas. Theoretical modeling results indicate that large 
view zenith angles (VZA) can lead to underestimation of fractional snow cover (FSC) by reducing the amount 
of the ground surface that is viewable through forest canopies, and by increasing uncertainties during the 
gridding of MODIS data. At the end of the MODIS scan line, the total modeled error can be as much as 50% 
for FSC Empirical analysis of MODIS/Terra snow products in four forest sites shows high fluctuation in FSC 
estimates on consecutive days. In addition, the normalized difference snow index (NDSI) values, which are 
the primary input to the MODIS snow mapping algorithms, decrease as VZA increases at the site level. At 
the pixel level, NDSI values have higher variances, and are correlated with the normalized difference vegeta- 
tion index (NDVI) in snow covered forests. These findings are consistent with our modeled results, and imply 
that consideration of view angle effects could improve MODIS snow monitoring in forested areas. 

© 201 1 Elsevier Inc. All rights reserved. 


1. Introduction 

Snow covers up to 50 million km 2 (over 30%) of the Earth's land 
surface seasonally (Dozier, 1989; Vikhamar & Solberg, 2003a). The as- 
sociated phase changes and unique high albedo of snow dominate 
local surface temperatures, impact radiative transfer and energy bal- 
ance, and influence global climate. Accurate snow extent mapping is 
valuable for various studies, such as surface runoff and water assess- 
ment, energy balance models, and meteorological models (Konig et 
al., 2001). 

Satellite sensors provide an efficient way to monitor snow as it en- 
ables mapping of large areas with high temporal frequency. Binary 
snow maps, which flag pixels as either snow-free or snow-covered, 
have been provided since 1966 when the first operational snow 
cover mapping was done by NOAA (Matson et al., 1986). Data from 
the MODIS instruments onboard Terra (launched in 1999) and Aqua 
(launched in 2002) are used to produce daily global snow cover prod- 
ucts at 500 m resolution (Hall et al., 2002). Some studies extended 
their focus to estimate the fractional snow cover (FSC), primarily 
through spectral mixture analysis (Painter et al., 2009; Romanov et 
al., 2003; Rosenthal & Dozier, 1996), a variety of models 
(Metsamaki et al., 2005; Vikhamar & Solberg, 2003b), and the use of 
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normalized difference snow index (NDSI) (Barton et al., 2001; 
Salomonson & Appel, 2004, 2006). 

Forests make remote sensing of snow more challenging especially 
for the MODIS sensor with wide scan angles (e.g. Hall et al., 2001; 
Painter et al., 2009; Salomonson & Appel, 2004,). The importance of 
accurate snow monitoring in forested regions is highlighted by the 
fact that forests cover approximately 40% of the seasonally snow- 
covered areas in North America (Klein et al., 1998). The presence of 
trees obviously complicates monitoring of snow using remote 
sensing as trees obscure snow on the ground surface. This effect is 
exacerbated by the wide view angles of MODIS as the viewable gap 
fraction through forests decreases as view zenith angles increase 
(Liu et al., 2008). In addition, MODIS gridded products have an 
inherent level of noise (Tan et al., 2006) that makes it difficult to 
provide accurate snow cover maps. 

There are other approaches for mapping the snow cover in forest- 
ed areas. Vikhamar and Solberg (2003b) in Norway used a linear 
spectral mixture model, which includes endmembers for snow, coni- 
fer, braches of leafless deciduous trees and snow-free ground, to cal- 
culate FSC. The validation of results from MODIS data using Landsat 
Enhanced Thematic Mapper Plus (ETM + ) images acquired the 
same days as the MODIS images showed good results when all forest- 
ed and non-forested areas are included. Given the use of Landsat for 
validation, the MODIS results are near nadir, which is the simplest 
scenario. It would be interesting to test this approach at high VZA. 
Metsamaki et al. (2005) developed a semi-empirical model with 
three contributors (wet snow, forest canopy and snow- free ground) 
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to estimate FSC from the reflectance observation. Their results for the 
Advanced Very High Resolution Radiometer (AVHRR) data are vali- 
dated against in situ measurements from weather stations and snow 
courses located in Finland. 

The purpose of this paper is to improve the understanding of 
MODIS snow mapping in forested areas and the MODIS snow prod- 
ucts. Specifically, we explored the effects of the wide view angles of 
MODIS with both model simulation and empirical analysis of 
MODIS/Terra products. The modeling process consists of two major 
steps: i) the Geometric-Optical Radiative Transfer (GORT) model is 
used to estimate the viewable gap fractions (VGF) through canopies, 
and ii) simulation of the gridding process of MODIS data (i.e. the pro- 
cess of taking MODIS swath observations and assigning them to the 
grid cells used for the MODIS products, see Wolfe et al., 1998 for de- 
tails). Empirical analysis of MODIS/Terra products is done to assess 
the view angle effects for four forested sites. 

2. Background 

2.1. Snow optical properties 

The underlying idea of using optical sensors to map snow is that 
snow has a distinctive spectral signature. The reflectance is very 
high in the visible wavelengths, drops steeply in the Near-IR and re- 
mains low for longer wavelengths (Wiscombe & Warren, 1980). 
Thus, NDSI, the normalized difference between green and Mid- 
Infrared reflectance, is widely used for characterizing snow (e.g. 
Dozier, 1989; Hall et al., 1995). 

Green -MIR _ MODIS Band 4-MODIS Band 6 
— Green + MIR — MODIS Band 4 + MODIS Band 6 ( ' ) 

Snow albedo is influenced by many factors, like grain size, finite 
depth, and impurities (Warren, 1982). Generally, grain size will in- 
crease with age and snow reflectance will decrease in all wave- 
lengths, but more at short wavelengths. Snow reflectance can vary 
from 0.9 for fresh snow to 0.5 or less for older snow in the visible 
wavelengths. Other factors that influence the reflectance of snow, 
like illumination angle, snow density, and liquid water content, can 
be considered equal in magnitude to the effect of changing effective 
grain sizes (Warren, 1982). Finite depth and impure mixtures also 
have considerable influence on snow reflectance in the visible wave- 
lengths, but little effect at 0.9 pm and beyond, where ice becomes a 
strong absorber (Wiscombe & Warren, 1980). 

2.2. The MODIS instrument and snow mapping algorithms 

The MODIS instruments onboard the Terra and Aqua satellites are 
multispectral imaging radiometers with 36 spectral bands from ap- 
proximately 0.4 to 14.0 pm (Barnes et al., 1998). The spatial resolu- 
tion of the MODIS instruments varies with spectral band, and ranges 
from 250 m to 1 km at nadir. MODIS provides daily, global coverage 
with a 2330 km swath width and a =b 55° across-track scan angle. 
The view zenith angle approaches 65° at the end of a scan line. 

The initial algorithm of MODIS binary snow mask used two criteria 
to identify snow pixels (Hall et al., 1995): NDSI greater than 0.40, and 
the reflectance in MODIS Band 2 (841-876 nm, near-infrared band) 
greater than 0.11. However, errors were found in forested areas 
where the vegetation masks snow on the ground. Therefore, Klein 
et al. (1998) improved the algorithm with the help of the Normalized 
Difference Vegetation Index (NDVI). In the new algorithm, individual 
pixels with NDSI values lower than 0.40 are identified as snow if NDVI 
values indicate that the pixels are forested. In addition, MODIS Band 4 
(545-565 nm, green band) should have a minimum reflectance of 0.1 
to prevent incorrect mapping of dark forests as snow. 


Salomonson and Appel (2004, 2006) developed the algorithm for 
FSC in Version 5 MODIS snow products (MO/YD10A1 ). It uses a linear 
relationship derived via regression to retrieve FSC from NDSI. In this 
operational algorithm, higher NDSI values lead to higher estimates 
of FSC. For Terra MODIS, the function is: 

Fractional Snow Cover = —0.01 + 1.45 x NDSI (2.2) 

This regression had high correlation coefficients and acceptable 
root mean square errors (RMSE) compared to FSC maps derived by 
aggregating Landsat TM binary snow maps. While this attempt to 
map fractional snow cover is validated with Landsat TM viewed 
from nadir, MODIS often collects data from much larger VZA and 
can lead to bias in heavily forested areas. 

2.3. View geometry effects on snow mapping in forested areas for MODIS 

View geometry may be one of the major error sources in MODIS 
snow mapping algorithms in forested areas. First, the gaps in forest 
canopies, which are essentially the detectable snow fraction in win- 
ter, will be lower at off-nadir views (Liu et al., 2008). Topography 
can also contribute to the view geometry effect, but is not discussed 
in this study. Second, toward the edges of MODIS swaths, the GIFOV 
(Ground-projected Instantaneous Field of View) become much larger 
than the grid cells to which they are assigned, with the results being 
that individual observations may cover several adjacent grid cells and 
significantly overlap each other (Wolfe et al., 1998). As a result, there 
are mismatches between swath observations and predefined grid 
cells, especially for observations at high view angles (Tan et al., 2006). 

3. Methods 

3.1. Modeling the reflectance of forests with snow 

The scenarios for our models include green trees and snow pre- 
sent under the trees. With bright clear atmospheres, snow rapidly 
melts off of boughs and fall below the canopy. We acknowledge 
there are times when snow remains on the boughs and can stay 
there for several days or weeks. However, we have not included 
these cases in our modeling as in these cases it will be much easier 
to find areas as snow-covered. 

The viewable gap fraction (VGF) at nadir equals the maximum 
viewable snow fractions (Liu et al., 2008): 

Max Viewable Fractional Snow Cover = VGF (Nadir) = 1— Forest Cover 

(3.1) 

The analytical form of the hybrid Geometric-Optical Radiative 
Transfer (GORT) model (Ni et al., 1999) was used to model the bidi- 
rectional reflectance of a forest stand with snow on the ground. In 
this model, the signal received by a sensor is a linear combination of 
the spectral signatures of four elements and their areal proportions 
(Li & Strahler, 1986). 

R = K c C + K g G + K t T + K z Z (3.2) 

where R is the total stand reflectance, K c ,K g ,K t ,K z are the areal propor- 
tions of sunlit crown, sunlit background, shadow crown and shadow 
background, respectively; and C, G, T, Z are the reflectances of these 
four components. 

This model combines the principles of geometric optics and radia- 
tive transfer theory. It treats the tree crowns as a collection of discrete 
ellipsoids that cast shadows on themselves and the background. The 
crown centers of these discrete ellipsoids are distributed randomly 
in space and uniformly between upper and lower height boundaries. 
The geometric optical method is applied to model the irradiance that 
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passes through between-crown gaps directly to the ground surface, is 
scattered by background, and exits without further scattering by veg- 
etation. Radiative transfer is used to calculate the probability of scat- 
tering based on path length through the canopy and provides the 
radiance that passes though and is scattered by plant crowns. 

Required inputs of the GORT model include spectral properties of 
leaves and snow, canopy parameters, and illumination angles. Table 1 
lists the typical values of spectral properties of leaves from the PROS- 
PECT model (Jacquemoud et al., 1996) and snow reflectance from the 
model of Wiscombe and Warren (1980). These values of reflectance 
were chosen to simplify this simulation step and to illustrate the ef- 
fects of view angles on the reflectance of canopies. Canopy parame- 
ters were from field measurements of the stands at the Fool Creek 
ISA and the St. Louis Creek ISA, which are both parts of the Fraser 
Meso-scale Study Area (MSA) in Colorado (Table 2). A more complete 
description of the canopy data can be found in Liu et al. (2004). Solar 
zenith angle was fixed at 50° for exploratory purposes, but obviously 
varies over time. 


3.2. Simulation of the MODIS gridding process 

The term “gridding” denotes the process used to allocate image 
observations into an output image grid (Wolfe et al., 1998). The 
MODIS gridding process puts MODIS swath observations into globally 
pre-defined grid cells. This process can introduce the “pixel shift” ef- 
fect, as the observations from any particular overpass and the grid 
cells are often misaligned (Tan et al., 2006). 

The MODIS gridding process was simulated to quantify the effect 
of view geometry on snow mapping accuracy. The ground reference 
data are a 30 m fractional forest cover map (Fig. 1) obtained from 
the USDA Forest Service Region 5— Remote Sensing Lab. This map 
covers an area of 147 x 147 km in the Sierra Nevada Mountains of Cal- 
ifornia, and classifies the forests into nine classes with 10% cover in- 
tervals (Table 3). Pixels with the value zero could either be pixels 
without data or pixels with 0% forest cover, thus they were excluded 
in further statistical analysis. 

To focus on the view angle effects, we assumed the landscape to 
consist of a single tree species and that all the forests have the same 
crown shapes as the St. Louis Creek stand. The only parameter that 
varies between pixels is the tree densities. Also we assumed the 
snow reflectance to be typical and the topography or atmospheric ef- 
fects to be minor. 

The method of Tan et al. (2006) was chosen to simulate the MODIS 
gridding process, and it took the following steps: 1 ) the VGF map was 
derived from the original 30 m forest cover map; 2) the scan angle, 
view angle, and MODIS GIFOV were calculated for each pixel with ar- 
tificial satellite positions; 3) we put the 30 m VGF map under the de- 
sired view angle, and all the 30 m pixels falling into individual MODIS 
observations were weighted by the MODIS point spread function to 
estimate the MODIS 500 m swath data; and 4) the MODIS 500 m 
swath data were placed into predefined grid cells based on the obser- 
vation coverage (Yang & Wolfe, 2001) to obtain estimates of the 
MODIS 500 m gridded VGF product. The results of our estimates 
were validated against the reference map at 500 m resolution 
(resampled from the original 30 m version). 


Table 1 

Spectral properties of leaves and snow for the MODIS bands. 



B1 

(645.5 nm) 

B2 

(856.5 nm) 

B4 

(553.6 nm) 

B6 

(1627.1 nm) 

Leaf reflectance 

0.05 

0.49 

0.10 

0.36 

Leaf transmittance 

0.03 

0.42 

0.07 

0.35 

200 pm Pure Snow 
Reflectance 

0.96 

0.85 

0.97 

0.04 


Table 2 

Canopy parameters as inputs in GORT model. 


Stands 

hi 

(m) 

h2 

(m) 

A 

(trees/m 2 ) 

R 

(m) 

b/R 

FAVD 

(m- 1 ) 

Dense stand in Fool Creek 

9.72 

21.97 

0.086 

1.53 

3.78 

0.578 

Intermediate stand in Fool 
Creek 

4.22 

14.07 

0.125 

1.13 

3.46 

0.495 

Sparse stand in Fool Creek 

3.45 

11.49 

0.116 

0.95 

3.31 

0.570 

St. Louis Creek 

7.75 

17.24 

0.058 

1.72 

3.09 

0.373 


hi : the lower boundary of crown centers; h2: the upper boundary of crown centers. 

crown centers density; R: horizontal crown radius; b: vertical crown radius, 
FAVD: foliage area volume density. 


In addition, we ran the GORT model with the calculated view ge- 
ometry information to add the view angle effects of forest masking, 
and then simulated the MODIS gridding process. This step helps to as- 
sess the total view geometry effects from both forest masking and 
MODIS gridding. 

3.3. Empirical analysis of MODIS /Terra products 

Four forested areas in the Northern US were picked for our empir- 
ical study of MODIS data from Jan 1st to May 19th in 2004 (140 days). 
Table 4 lists the precise location and meteorology information 
(WWW1). Since all these sites have frequent snowfall and freezing 
temperatures that last for the first 90 days, we assume that these 
areas have continuous snow cover during this time period. Fig. 2 
shows the four areas at high resolution from Google Earth™, includ- 
ing: 1) a large flat region in Eagle River, Wisconsin; 2) a relatively 
flat area in Springfield, Vermont; 3) Harvard forest in Massachusetts; 
and 4) a small mountain area close to Lebanon, New Hampshire. All of 
these sites are almost homogeneous and fully covered by forests. For 
each location, we chose a 3x3 pixel (1.5 x 1.5 km) window rather 
than single pixel to minimize the influence of pixel-to-pixel 
variability. 

The FSC data were obtained from the MODIS/Terra snow cover L3 
daily global product at 500 m resolution (MOD10A1). The MODIS/ 
Terra 500 m daily surface reflectance product (MOD09GA) was used 
to study the relationship between NDSI and sensor angles. NDSI 
values for each window were calculated from the averaged MODIS 
Band 4 (green) and Band 6 (mid-infrared) reflectances at 500 m res- 
olution. Because snow-covered forests usually have a much lower re- 
flectance in the green band than pure snow, we set a minimum 
reflectance threshold of 0.01 in Band 4 to eliminate unrealistically 
low reflectances, which could result in very low NDSI values. 
MOD09GA is delivered with 1 km sensor and solar angle data and 
quality data, which were resampled to 500 m resolution. 

The quality data provide internal cloud and cloud shadow flags, 
but only the internal cloud flag was used as the MODIS cloud mask 
is rather conservative (Stroeve et al., 2006). In addition, if there was 
only one clear observation in the 3 x 3 window, the reflectance data 
for that day would not be used. The results below show that some 
pixels could not be captured by the internal cloud algorithm, especial- 
ly when they were mixed or spectrally similar to the clouds. Since 
even thin clouds can influence the NDSI strongly, we visually checked 
all the available days after applying the internal cloud mask to make 
sure they are free of clouds. 

4. Results 

4.1. Modeled view angle effects of canopy masking 

VGF decreases consistently with decreasing VZA for all four sites 
(Fig. 3). The viewable background of all four stands at nadir drops 
by about 50% at the largest MODIS view zenith angle. For instance, 
the Dense Stand has a 60% gap at nadir but only less than 4% gap at 
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Fig. 1. The 30 m fractional forest cover map of the Sierra Nevada Mountains in California. 


VZA of 65° (nominally, MODIS has ±55° scan angle, which is about 
±65° VZA). 

We chose the St. Louis Creek stand as an example to show GORT 
modeled reflectances. Reflectance in the principal plane is influenced 
by the hotspot effect, where the view geometry coincides with the 
solar geometry. Considering large solar zenith angle in winter, 
MODIS sensor rarely collects data near the hotspot, so results in the 
cross-principal plane are closer to the real MODIS acquisitions. In 
Fig. 4, modeled reflectances generally decrease in the visible bands, 
but increase in the mid-infrared and near-infrared as view angles in- 
crease. For landscapes with only snow, mid-infrared reflectance 
should be much lower than green; but here we see that the trees 
have a significant effect: mid-infrared reflectance is similar to green 
wavelengths and even greater at large VZA. 

The dependence of NDSI on the VZA (Fig. 5) illustrates how the 
view angle can influence the performance of the current MODIS 
snow algorithm for snow-covered forests. The NDSI and NDVI 
values were calculated from GORT modeled band reflectances. These 
values of all four stands have similar dependence on the VZA, but 
with different magnitudes. Although the NDSI values remain 
positive for the Sparse Stand, NDSI decreases rapidly with 
increasing VZA, while the NDVI values increase slightly. The MODIS 
snow mask can correctly identify these stands as snow-covered if 
viewed at nadir, but would have troubles at large VZA. Snow cover 
in the Dense Stand will be missed when the VZA is greater than 5° 
because the NDSI values are negative. NDVI remains positive at all 
view angles, even for the Sparse Stand at nadir. Since the NDSI 
values change with VZA, the MODIS snow products may give 
different fractional snow cover values for observations from 


Table 3 

Class information of the forest cover map in Sierra Nevada, California. 

Forest cover (%) 0 10 20 30 40 50 60 70 80 90 

Percentage of all 48.2 1.2 2.6 4.9 7.4 7.6 6.3 5.3 5.7 10.6 

pixels (%) 


different days, even though the snow cover does not change 
significantly. 

The important factors in the GORT simulation include the shape of 
the trees, forest and snow optical properties, and solar geometries. 
The ellipsoid shape was used in this study instead of a conical 
shape, since it has already been shown to be effective for estimating 
VGF though forest canopies (Liu et al., 2008). In addition, we modeled 
the reflectance of individual bands and NDSI against VZA with varying 
tree and snow optical parameters at different solar zenith angles. 
Spectral reflectance shows strong dependence on some of these in- 
puts, for example, leaf transmittance, snow grain sizes and impurities, 
and solar zenith angles. However, the outputs of our simulations still 
show the strong correlations between NDSI and VZA presented here. 

42. The view angle effect of the MODIS gridding process 

The symmetric pattern in Fig. 6 indicates that the VZA introduces 
random uncertainties during the MODIS gridding processes because 
the GIFOV of MODIS sensor is changing with increasing VZA. This ef- 
fect could result in either overestimation or underestimation of the 
snow cover fraction, depending on the forest cover in the areas sur- 
rounding the grid cells. VGF estimates have higher uncertainty from 
the MODIS gridding process at larger VZA, with RMSEs of 6.46%, 
8.15% and 11.1% for a VZA of 5°, 30°, and 55°, respectively. 

Table 4 

Information about four forest sites. 

Location Coordinates Jan 1st 2004 to Mar 30th 

2004 (first 90 days) 

Total snowy Avg max 


days temperature 

Eagle River, Wisconsin 45°55'39"N 88°50'14"W 61 -2.8 °C 

Springfield, Vermont 43°23'05"N 72°32'18"W 40 -0.56°C 

Harvard forest, 42°31 / 04"N 72°12 / 49"W 31 0.56 °C 

Massachusetts 

Lebanon, New Hampshire 43°43'06"N 72°08'42"W 41 -1.1°C 
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Eagle River, Wisconsin Springfield, Vermont 



Harvard forest, Massachusetts Lebanon, New Hampshire 



Fig. 2. A visual view of four study sites from Google Earth™. The white squares in maps are approximately 1.5 x 1.5 km, which represent the 3x3 pixel windows in MODIS 500 m 
gridded products. 


The effects of combing both changes in VGF and the gridding pro- 
cess are shown in Fig. 7. Increasing VZA can substantially reduce over- 
all VGF and gradually shifts the pattern to underestimation. The 
symmetrical pattern at nadir indicates that the effects are mainly a re- 
sult of the gridding process. At off-nadir view angles, the canopy 
masking effect dominates. In Fig. 7c, most VGF estimates are lower 
than 10%, which is the lowest value in the original map. The RMSE 
values also increase with increasing VZA, and could be as high as 
31% at VZA of 50°, or close to the end of MODIS scan lines. 



Fig. 3. Viewable gap fractions (VGF) as a function of view zenith angle (VZA) for the 
four stands modeled with the GORT model. The black solid line on the right represents 
the end of the MODIS scan line. 


4.3. Empirical analysis 

When the averaged NDSI values for the four forested sites are 
plotted as functions of view zenith angle, the effect of high view an- 
gles in forested landscapes is obvious (Fig. 8). Second-order polyno- 
mial curves can fit these points well. As observations move away 


St. Louis Creek 



Fig. 4. GORT modeled band reflectance as a function of view zenith angle in the cross- 
principal plane for the St. Louis Creek stand. Solar zenith angle is fixed at 50°. 



Q, Xin et al. / Remote Sensing of Environment 118 (2012) 50-59 


55 



Fig. 5. GORT modeled NDSI and NDVI as a function of view zenith angle in the cross-principal plane for four different stands. Solar zenith angle is fixed at 50°. 


from nadir, NDSI decreases at all sites. Because the presence of forest 
can lower the green reflectance significantly, NDSI can be negative 
even with snow on the ground. For example, at the Springfield site, 
NDSI ranges from about 0.5 at nadir to -0.3 off-nadir. After 
snowmelt, NDSI remains low and is independent of view angles. 
These findings are consistent with the GORT model results and 
reflect the view angle effects of canopy masking (Fig. 5). In addition, 
the different rates of decline may reflect canopy densities. The 
polynomial curves at the Eagle River, Springfield, and Lebanon sites 
are close to the modeled curves of the Intermediate stand, while the 
Harvard forest site is similar to the Dense Stand. 

Application of the MODIS snow mapping algorithm needs to be 
done at the individual pixel level, rather than at the aggregated 
level shown in Fig. 8. The NDSI values of individual pixels (Fig. 9) 
have more variability than the averaged NDSI of the window at the 
same site, since observations from a single pixel usually couple with 
the MODIS gridding effect. This could explain why some observations 
are far from the fitted lines for the 3 x 3 window. After snowmelt, the 
NDSI values of both the center pixel and the 3x3 window are consis- 
tently low. The MODIS snow mask performs fine for VZA up to about 
30° at the Eagle River, Springfield and Lebanon sites, but is ineffective 
at the dense Harvard Forest site. Applying the MODIS fractional snow 
regressions, the FSC could range from zero to around 80% even if all 
these sites remained snow-covered during the snow season. 

Fig. 10 shows the FSC values in the MODIS snow cover products for 
the four forest sites. The dramatic changes of the FSC values on 


consecutive days are unrealistic during the snow season. For example, 
the FSC decreases from over 40% to about 0% and then increases to 
50% in just 7 days at the Springfield Site. The FSC estimates are always 
zero after snowmelt (near Day 100) and can easily be separated from 
the snow season. These findings cannot be explained by the snow- 
melt or snowfall events, but are consistent with our expectations for 
significant view geometry effects: MODIS has very different view an- 
gles for the same area on consecutive days, thus the fluctuating FSC 
estimates are results of both underestimation of canopy masking ef- 
fects and uncertainties in the gridding process. 


4.4. Correlation between NDSI and NDVI 

Strong correlations are found in the GORT model results when 
NDSI are plotted against NDVI with varying view angles in the 
cross-principal plane (Fig. 11). Empirical analysis also shows that 
NDSI is highly correlated with NDVI in snow-covered forests 
(Fig. 12). Linear regressions between NDSI and NDVI are performed 
separately for each site with all the cloud-free data of each pixel in 
the 3x3 window, excluding those after snowmelt. The four sites 
show high R-square values of 0.72, 0.80, 0.80, and 0.75, respectively. 
This is another way of showing the strong view angle effects in 
MODIS: when viewing away from nadir, we see more trees and less 
snow, thus NDSI deceases while NDVI increases. This finding also im- 
plies that NDSI is probably not a good index for deriving FSC from 
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Fig. 6. Comparison of VGF ground reference with estimation. The model estimations come from the simulation of the MODIS gridding processes. Colors represent the numbers of 
pixels. The solid line is the 1:1 line, where the estimates are equal to ground references. 
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Fig. 7. Same as Fig. 6, except that model estimations are simulated by both the GORT model to include the forest masking effect and the MODIS gridding process. Notice that as the 
view zenith angle increase, so do the magnitudes of errors. 


MODIS products for forested areas, since it is so closely related to 
NDVI. 

The percentages of observations identified as snow-covered by the 
MODIS binary snow algorithm are shown in Fig. 12. The observations 
that fall in the triangles or to the right of the vertical line (NDSI> 0.4) 
in Fig. 12 will be classified as snow in the latest version of the MODIS 
binary mask, whereas only pixels to the right of the vertical line 
would have been classified as snow in the old version of the snow 
mask. The initial binary snow algorithm cannot effectively detect 
snow in forested areas. The added triangle in the NDSI-NDVI 
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domain improves the snow mapping accuracy by recognizing that 
additional 28%, 25%, 4% and 32% of these stands are snow-covered. 
However, over 60% of the observations are still not correctly 
identified as snow-covered, especially those acquired at large VZA. 
Pixels after snowmelt are generally very low in NDSI and 
distinguishable from the others, although some snow-covered pixels 
may also exhibit very low NDSI values. While NDVI stays above zero 
for forested areas, NDSI is often negative even when snow is present. 
This implies that the usage of NDSI alone to discriminate snow is 
quite limited in snow-covered forests. 


Springfield 
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Fig. 8. Site-level NDSI values as a function of view zenith angle for four selected forest sites. NDSI values are calculated from averaged band reflectance in a 3 x3 pixel (1.5 x 1.5 km) 
window. The solid line is a second-order polynomial curve fitting the values during the snow season. 
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Fig. 9. Pixel-level NDSI values as a function of view zenith angle for the four selected forest sites. NDSI values are calculated from the band reflectance of the center pixel in a 3 x 3 
pixel (1.5 x 1.5 km) window. The solid line is a second-order polynomial curve fitting the site-level values during the snow season. Dots represent the days with snow that are de- 
tectable by the MODIS binary snow algorithm, circles represent the days with snow that cannot be detected by MODIS binary snow algorithm, and asterisks represent the days after 
snowmelt. 


5. Discussion and conclusions 

We have demonstrated why it is difficult to monitor snow accu- 
rately in forested areas for sensors with wide view angles like 
MODIS. Both model outputs and empirical research show that the 
snow cover identification and fractional snow coverage using the 
MODIS algorithms vary considerably in forested areas. The variation 



Fig. 10. Fractional snow cover from MODIS snow products (MOD10A1) at four forest 
sites from Jan 1st to May 19th, 2004. The FSC values of each site are from the center 
pixel in a 3 x3 pixel (1.5 x 1.5 km) window. 



Fig. 11. GORT modeled NDVI against NDSI for the four forest stands with snow on the 
ground. NDVI and NDSI values are calculated from reflectances in the cross-principal 
plane. For each stand, nadir observations are at the end of the line at the bottom 
right. As view zenith angle increases, the values move to the upper left with a maxi- 
mum view zenith angle of 65°. The MODIS snow mask flags observations as snow if 
they fall in the triangles or to the right of the vertical line (NDSI >0.4). 
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Fig. 12. NDVI against NDSI for the four forest sites. Dots represent observations with snow cover, while circles denote the days after snowmelt. The MODIS snow mask flags obser- 
vations as snow if they fall in the triangles or to the right of the vertical line (NDSI > 0.4). The percentages illustrate the proportion of observations that would be classified as snow/ 
no snow in the MODIS snow mask. The functions are the linear regressions between NDVI and NDSI in snow season, along with key statistics. 


in the MODIS snow products is largely due to the widely varying view 
angles from the MODIS instrument on a day-to-day basis. This study 
helps us better understand the quality of MODIS daily snow cover 
products, and indicates that consideration of the VZA is required in 
forested areas. 

The modeling process showed that observations from large VZA 
can result in significant underestimation of VGF through forest cano- 
pies, and consequently lower the value of NDSI for winter forests. Fur- 
thermore, the MODIS gridding effects will introduce uncertainties in 
observations at large scan angles. When planning future satellite sys- 
tems for operational snow monitoring, it should be noted that in- 
creased frequency of near nadir observations will improve results. 

Empirical analysis of the MODIS products showed high day-to-day 
variability of the FSC for forested areas during a period with consis- 
tent snow cover. The daily averaged NDSI for a 3x3 window of 
snow-covered pixels at four forested sites were found to decline 
with increasing VZA (with some NDSI values being negative). The sin- 
gle center pixels of the 3x3 windows exhibited higher variability 
than the windows at the same site due to uncertainties introduced 
by the gridding process. There were strong correlations between 
NDSI and NDVI for four study sites. These findings support the 
model results, and confirm the implications of the view angle effects 
on viewable snow fractions. Further, they imply that it is difficult to 
derive snow fractions solely from NDSI for snow-covered forests, be- 
cause vegetation obscures snow and considerably lowers NDSI. 

This research indicates the potential for improving snow mapping 
in forests by explicitly taking into account the view angle effect. 


Although the MODIS binary snow algorithm improved its perfor- 
mance in forested areas by integrating NDVI, it still has omission er- 
rors for observations at large VZA or in dense forests. In addition, 
the MODIS fractional snow products have trouble in estimating the 
FSC in forests, since NDSI is strongly view angle dependent and close- 
ly related to NDVI. It also appears that adding NDVI into consideration 
would help the FSC algorithm for the MODIS sensors. 

A flexible way to compensate the view angle effect is to derive the 
function between VZA and NDSI or FSC in forested areas. This is not an 
easy task because that relationship is also dependent on the forest 
cover and the canopy structure. Dozier and Frew (2009) provided 
an alternative method to interpolate the estimates of the daily snow 
cover, thus smoothing the effects from sensor viewing geometry. 
The simplest solution may be to limit the use of MODIS snow data 
to VZA in the range of 0° to 30° if pixels can be identified as forests. 
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